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Abstract 

Hydrodynamics predicts long-lived sound and shear waves. Thermal fluctuations in these waves 
can lead to the diffusion of momentum density, contributing to the shear viscosity and other 
transport coefficients. Within viscous hydrodynamics in 3+1 dimensions, this leads to a positive 
contribution to the shear viscosity, which is finite but inversely proportional to the microscopic 
shear viscosity. Therefore the effective infrared viscosity is bounded from below. The contribution 
to the second-order transport coefficient r,r is divergent, which means that second-order relativistic 
viscous hydrodynamics is inconsistent below some frequency scale. We estimate the importance of 
each effect for the Quark-Gluon Plasma, finding them to be minor if ij/s = 0.16 but important if 
rj/s = 0.08. 
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I. INTRODUCTION 



Heavy ion collisions at RHIC [1-4] and the LHC [5, 6] produce a medium whose evolution, 
at least at early times, is well described by hydrodynamics with a very small viscosity [7-15] 
(at least when normalized to the entropy density [16]). A major theoretical goal is now to 
determine this viscosity as accurately as possible, by modeling the development of a heavy 
ion collision with viscous hydrodynamics. In order to do so it is necessary to go beyond the 
"first-order" (Navier-Stokes) formalism, because the first-order formahsm leads to acausal 
and unstable evolution in a relativistic setting [17-21]. As shown by Israel and Stewart 
[22-24] , this problem can be cured by working instead with hydrodynamics expanded to the 
second order in gradients, which should have the added advantage of being more accurate. 
Several groups have been involved in studying the hydrodynamics of heavy ion collisions 
using such second-order formalisms [10, 25-30]. However, within the community which 
studies hydrodynamics and kinetics of atomic gases, it has been known for almost 40 years 
that the gradient expansion in hydrodynamics fails beyond the first order [31]. Is this 
also true in the relativistic setting? If so, what implications does it have for the study of 
hydrodynamics via the second-order formalism? 

II. SETUP AND INTUITIVE ARGUMENT 

We start by reviewing relativistic hydrodynamics to second order. Hydrodynamics is the 
modeling of a fluid by solving the stress-energy conservation equations^ 

d^T'^'^ix) = (2.1) 

assuming some functional form for the stress tensor. Ideal hydrodynamics assumes the 
equilibrium form, 

T^^-{x) = {e{x) + P) M^(a;) ^"(a;) + Pg'"' , P = P(e) (2.2) 

where Qfj^i, — Diag[— 1, -|-1, -|-1, -|-1] is the metric, is the 4-velocity determining the rest- 
frame (normalized so u^u^ = —1), e = u^u^T^^ is the rest-frame energy density, and P 
is the pressure as determined by the equation of state P = -P(e). Viscous hydrodynamics 
assumes that the fluid is near equilibrium so that T'^" is close to this form. Assuming that 
equilibration is a fast, local process, corrections to this form can be written in terms of an 
expansion in gradients. Israel and Stewart showed [23] that a slight re-organization of the 



If there are conserved currents one also considers current conservation d^Jj^ = 0. In equilibrium 
Ja = n-aU^, out of equilibrium there can be derivative corrections. However, in ultra-relativistic heavy 
ion collisions the density of the conserved baryon number in the central rapidity region is small, so we 
will neglect it and will not discuss conserved currents further. 
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second-order derivative expansion yields stable equations which are correct to second order 
provided a particular additional term appears at second order; 



T^^ix) = + W , W = -2ri,id^''u''^ + (u'^dall'"' + in^^a^w") (+ other terms) 

^ (2.3) 

where rjci is the "classical" viscosity coefficient that one would obtain from a microscopic 
calculation using the Kubo formula based on (2.3). The angular brackets in d^^u"'^ mean 
that the indices are to be symmetrized, projected to be spatial in the frame given by 
and trace-subtracted.^ The extra terms include bulk viscosity and nonlinear effects, and are 
catalogued in [32-34]. To simphfy the discussion here we will consider a conformal fluid, in 
which case the bulk viscous term is absent. While this is not a very good approximation 
for QCD near the transition temperature, in practice the inclusion of bulk viscosity and a 
realistic equation of state would have only a small influence on our final results.^ 

Solving the hydrodynamic equations for small fluctuations inu^, e about constant values, 
one finds two sorts of long-lived wave solutions, sound waves and shear waves. At lowest 
order in k,uj <^ ^/Vch they obey dispersion relations of 

Each sort of wave decays with time, but with a decay rate which vanishes quadratically in 
the small k limit, as well as becoming small for small rjci- Generically, out of equilibrium such 
waves will be present with large amplitudes. But equipartition of energy says that even in 
equilibrium such waves will be present, carrying energy which averages to T/2 per degree of 
freedom. The long decay times of these waves contradict the assumption that all degrees of 
freedom in a fluid equilibrate via rapid local processes. This imperils the assumption behind 
the gradient expansion in Eq. (2.3). If there are arbitrarily slowly equilibrating degrees of 
freedom, Eq. (2.3) can contain terms nonanalytic in gradients. In the nonrelativistic setting 
it is known that precisely this happens [31]. 

To see more intuitively how hydrodynamic waves can contribute to hydrodynamic co- 
efficients, consider a shear wave, as illustrated in Figure 1. Hydrodynamics says that this 
wave configuration will decay with time as the x-momentum carried by the fiuid diffuses 
in the y direction into neighboring regions, which are not flowing. The rate of this diffu- 
sion is controlled by shear viscosity, which is by deflnition the diffusion coefficient for the 
component of momentum transverse to the diffusion (the y-diffusion of x-momentum in this 
illustration) . But one of the mechanisms which can transmit x-momentum is hydrodynamic 
waves, such as sound waves, with wave lengths shorter than the hydrodynamic structure 



^ That is, defining P^" = g^^ + u^v^ which is a projector to local rest-frame spatial components, 2d^'^u"^ = 

^ In particular our results are dominated by the effects of shear waves, which are not sensitive to bulk 
viscosity or the equation of state. 
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FIG. 1: A shear wave; fluid moves to the left in a band of fluid near the middle of the figure. 
Viscosity determines the loss (by diffusion) of the forward motion of this fluid. Sound waves 
(dotted, in red) leaving the left-moving fluid carry, on average, net left-moving momentum, which 
is not compensated by sound waves arriving in the left-moving fluid. Hence sound waves contribute 
to viscosity. 

considered. Sound waves leaving the a;-moving fluid carry net rc-momentum away, diffusing 
away this component of momentum. The phase space of such waves scales as d^k which 
is very UV dominated; but the distance propagated before dissipation scales as 1/k^ as we 
just saw, leading to a contribution to viscosity which scales as ~ (Pk/k"^. This is IR flnite, 
though it would not be in 2 spatial dimensions [35]. The contribution is larger for smaller 
T^ci, both because hydrodynamic waves propagate further and because the range of k where 
hydro is valid expands at smaller r]ch so the phase space of hydrodynamic waves is larger. 

Similarly, can be interpreted as a relaxation time; if a fluid suddenly develops shear 
flow, is the time scale for momentum diffusion to be established. We already saw that 
hydrodynamic modes contribute to momentum diffusion, with a mode of wave number k 
contributing of order 1/k'^. The time scale for this mode to leave equilibrium and establish 
its contribution to momentum diffusion also scales as l/k"^. This suggests that such a mode 
contributes to by an amount proportional to leading to an / d^k/k'^ contribution to 
x^-, which is small k divergent. 



III. COMPUTATION OF THE CONTRIBUTION FROM HYDRODYNAMIC 
WAVES 

Ref. [32] derives the following Kubo relation for the shear viscosity rj^x and the relaxation 
time r,r in terms of the retarded correlation function for two T^^ stress-tensor operators: 

G™(a;, K) = P- tcorj^i + {vaTn ' ^) ' + • (3-1) 
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Here k is another transport coefficient discussed in [32]. We see from the above that k 
can be extracted from the zero frequency behavior of the stress-stress correlation function, 
G^^^^{uj — 0,kz)- The zero frequency retarded function equals the Euclidean correlation 
function, so k is a thermodynamic property [36], which is not sensitive to long- wavelength 
hydrodynamic waves. Since we are not interested in the value of k, we will consider G^'^^ 
at vanishing external spatial momentum A; = and small nonzero frequency u. 

To compute the contribution of hydrodynamic waves to rj and r^r, we compute their 
contribution to the above correlation function, along the lines of what Kovtun and Yaffe 
did for the symmetrized correlator [37]. The operator T^^ to be used in the Green function 
above is the hydrodynamic one as described in Eq. (2.2) and Eq. (2.3) above, allowing for 
thermally occupied fluctuations in and e. In the absence of fluctuations — (1, 0, 0, 0), 
€ — €o, P — €o/3, where eo is the equihbrium energy density. Therefore arises at order 
(Su) due to the viscous terms, and at order (5it)^ and higher from the equilibrium and 
viscous terms: 



rpxy 



-r/ci(a^M^+a^M'^)] + [(e+P)M^M^-r/ci(M^aoM^+M^aoM^)-5e^(a^«J'+a^M^)]+C'((5=') 



(3.2) 

The contribution at first order in fiuctuations arises only from the viscous term —r]c\{d'^u'^ + 
d^u^). The symmeMzec? correlator is 

G^^^^[l-order] = j xdte^^-'^ '' {( - Vaid'^uy + a%^](x, i)), ( - r]^\[d^uy + a%^](0, 0))}^ , 

(3.3) 

which is related to the retarded correlator via Gs{t^) = ^[1 + 2n5(cj)]Im G-r{uj). The ex- 
pression above is automatically 0{k^) and so it does not contribute in the /c — )■ limit. 
But there can also be local contributions, that is, contributions proportional to (5^(x) or its 
derivatives, arising from the overlap of T operators [contact terms] . We will present a direct 
calculation of such contact terms within hydrodynamic theory in a future publication [38]. 
For our current purposes we will instead extract them by using stress conservation, 

^^T^'y^Q =^ klG''y''y{k^,uj)^uj^G^y^y{k^,uj) (3.4) 

together with the known expression for Gg^°^ (see Eq.(33d) of [37]): 

ar'it,,.) = 2,..T(,-^j§0pl-,). (3.5) 

The first term in (3.5) is the contact term, which gives rise to the —iurjci term in Eq. (3.1). 
The second term is the contribution from the correlator of the —rjcid^u^ part of T"^^, which 
vanishes at small k. 

Now we extend this calculation to second order in fiuctuations. At this order there are 
higher order corrections to the terms involving a single power of 5u in T^^; and the lowest 
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order contribution from terms involving two powers of u^, u^, 5e in T^^. But since the terms 
involving a single power of Su in T"^^ always involve spatial derivatives, these terms all vanish 
in the /c — ?> limit, like the second term in Eq. (3.5) above. So we skip their calculation 
and concentrate on contributions from second-order in fluctuation terms in T^^. There 
are several such terms in Eq. (3.2) but we will concentrate on the term (e+P)-u^-u^, and 
explain why the other terms can be neglected at the end. Since u^u^ is already quadratic 
in fluctuations, we can neglect fluctuations in e and replace (e+P) with cq+Po = f^o- Since 
fluctuations in e will play no further role in the discussion we will henceforth write e, P for 



where ()s,r indicate whether the operators are to be symmetrized (S) or if the retarded 
correlator is to be used (R). Using the expression wc have for the stress tensor above, 
the symmetrized correlation function of two stress tensors at vanishing external spatial 
momentum is 



Gri2-order](a;,A; = 0) ^i^ + Pfj^ J -0^^ {Ar{u',p)Af(u - u',-p) 

+ A^^(a;',p)Ar(a;-a;',-p)) , (3.7) 




where we assumed small, nearly linear hydrodynamic fluctuations and small frequencies 
(details will be given in Ref. [38]). 

We could use the KMS relation Gs{uj) = —[1 + 2n5(cj)]lm Gb,{uj) to extract the shear 
viscosity from this correlation function, but since t^^ depends on the real part of the retarded 
function we would need to invert the KMS condition through a Kramcrs-Kronig relation to 
get the real part of Gr. It is more economical to compute Gr directly; 



G'r^[2-order](a;,^ = 0) = + / ^ / (^ ^ 

(Ar (cu', p)A^ico-co', -p) + A-^(c.', p)Ar (a;-a;', -p) 

+ A-(^', p)Af (a;-a;', -p) + A^^(^', p)Ar (^-o;', -p)) (3.8) 



which in the smaH frequency limit reproduces Gg^^^ when using the KMS condition. This 
contribution is to be added to the flrst-order contribution, which, as we discussed, reproduces 
the terms present in Eq. (3.1). 

Now Ag (cj', p) was determined by Kovtun and Yaffe [37]. They included the effect of the 
viscous rici term (first order gradients) but dropped the term (second order gradients), a 
procedure we will follow. One finds 



eo,^(eo)- 

Define the correlation function of the fluid velocity to be 



Agi,(^',p) = I dtd'^ e-^P-^+^'^*(«^(t,x)«^(0,0))s,R, 



(3.6) 





e + P [ p2 (^2 _ p2/3)2 + (^^p2^) 



2T IriP^Oj"^ 



2 + (^"' - 



(3.9) 



with 



~ e + P ' ~ 3(e + P) ■ 
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The first and second terms in A*-^ represent sound and shear waves respectively. This 
expression differs shghtly from the one found by Kovtun and Yaffe [37] because they assume 
-C 1 which allows them to split the sound mode contribution into separate terms for 
cu ~ p/ and u ~ —p/ \/3 propagation. 

In the small frequency limit, the factor 2 T in the numerator of Eq. (3.9) should really be 
interpreted as cu(l + 2nb{uj)), with nb{u}) = [e^^'^ — 1)~^ as usual. Then the KMS relation 
A^(p, cj) = —(1 + 2nfe(c(;))Im A^(p, cu) and the analytic properties of Ar (no poles in the 
upper half-plane) uniquely establish 



e + P 



pY ^ f^ij _ pY \ -7^p2 



i%p'^uj + {00"^ — p'^/3) \ p^ ) —iu + ^jjp'^ 



(3.10) 



We are now ready to compute Eq. (3.8). Since our calculation is done entirely within 
hydrodynamics, or equivalently assuming momenta to be small, our result will not be ap- 
plicable at large momenta, neither in the argument of G^^^, nor for the hydrodynamic 
propagators A*-' inside the integral of Eq. (3.8). Hence, we need to restrict the calculation 
to the highest wave number Pmax) or the inverse of the shortest length scale, where the hy- 
drodynamic description is valid. Cutting off the p integration at Pmax and considering first 
the shear-shear contribution we can perform the index contractions and angular integrals, 
finding 

14 1 /-Pmax ^ ^ r du' -I-.jrl 



14 1 /"Pmax r 



27r {"friP^ — ioj') (7r,p^ -|- ioj') 7^p^ — iu + iu' 
30^ Jo - ioj/C^lv) 

r,,2/r4^2^ 1 

(3.11) 



'max 



7T fP^ , 







2 iu ^7(47') 
-p - h 



27^ p^ - iuj / 



The p^ term is an uninteresting contribution of hydrodynamic waves to the pressure.^ The 
two "interesting" terms are 

GZl^-.^M ^ -^^S^ + + 1)^^-^ • (3-12) 

607r^7^ 2407r7| 

Comparing with Eq. (3.1), we see that the first term is a positive contribution to the shear 
viscosity arising from shear waves. The imaginary part of the second term is a frequency 
dependent reduction in the shear viscosity, which vanishes as a; — )■ 0. Therefore there is no 
problem defining the shear viscosity in terms of the zero frequency limit of dG'^'^'' /du. The 
real part of the co^^"^ term has the same sign as the term in Eq. (3.1), but the wrong u) 
dependence. This term can be interpreted as a frequency dependent correction to which 



The pressure contribution has an unexpected sign. In a weakly coupled theory, the contributions of 
ordinary particles in loops also have the wrong sign, which is over-canceled by a contact term, see [39] . 
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diverges at small frequency; or it can be interpreted as a breakdown of the validity of the 
hydrodynamic expansion beyond one- derivative order. 

We should also include the sound {p^jP /p^) terms in A*-' shown in Eq. (3.9) and Eq. (3.10). 
If we assume PmaxTr? <^ 1 then the sound part in (3.10) may be approximated as 



+ (a;2 - p73) 2 \uj + i%/2p'^ - vl^ ^ + ^%/2p2 + vl^ ) 

which has the advantage that now the bj' poles in (3.8) are simple. The integration is then 
straightforward and for the mixed shear-sound term one has 

57r2 i-p^ ^^p2 _ 3^(^ _ p/V3)/57, 3712^"^-- + ' 371^ "^^-^ ^ " " " " 

One thus finds that the mixed shear-sound term is suppressed with respect to the result 
(3.12) by extra powers of Pmax7?j- The sound-sound term contains a part that has the same 
structure as the shear-shear contribution, as well as other parts that are again suppressed 
by powers of PmaxT??- Evaluation of the "interesting" contribution to leading order in PmaxT?? 
gives 

GT\^<^V^..<^%')^-^i^i^^{i^r)u^^^ +0(p^a.7^,^^) (3.13) 

This is our main result. 

Let us finally comment about whether or not we need to consider other contributions to 
Q^^y coming, for instance, from the —rjciu^dou''^ term in Eq. (3.2). Since hydrodynamics 
is predicated on the convergence of the derivative expansion introduced in Eqs. (2.2,2.3), 
Pmax should be chosen as the largest momentum scale where successive terms in the series 
are successively smaller, which requires Pmax <^ (e + P)/Vci (comparing the zero and one- 
derivative terms), Pmax ^ Vci/ (VciTtt) — (comparing the one-derivative and two-derivative 
terms). This also ensures that hydrodynamic waves with p < p^ax will have Imu <^ p, so that 
hydrodynamic waves are well-defined, long-lived excitations in the plasma. And it ensures 
that the real parts of the propagating frequencies of the two sound waves and the shear 
wave are more widely separated than their imaginary parts. These conditions ensure that 
contributions to G'^^^ from higher-derivative terms in Eq. (2.3), and contributions arising 
from interference between sound and shear waves in Eq. (3.8), are small compared to the 
terms we have computed. In particular, terms arising from T^^ D — r^ciW^^oM^ in Eq. (2.3) 
will give rise to an integrand similar to Eq. (3.11) but with an extra power of {p^rif'- The 
resulting term will be analytic in the frequency and will give corrections suppressed by Pmax7^ 
relative to the terms we have computed. We will not attempt to compute such suppressed 
corrections here. 
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IV. INTERPRETATION AND DISCUSSION 



A. Viscosity 

We found that, besides the "classical" viscosity rjci, there is an additional contribution 
to viscosity as measured on very long distance and time scales, generated by relatively 
short-wavelength hydrodynamic waves. It is given by the coefficient of —iou in Eq. (3.13): 

where higher power corrections in Pmax can be neglected as long as Pma^Jri <^ 1- The new 
contribution from Eq. (3.13) scales as an inverse power of t^cI- Therefore ~ t^cI + 0{l/rici) 
has a positive minimum at a nonzero value of f]ci- This places a lower bound on the total 
(infrared) value of rj in any given theory. 

To estimate the size of the new term and the minimum possible viscosity, we need to 
estimate Pmax- As discussed, pmax should the the largest wave number such that the gra- 
dient expansion converges and a hydrodynamic description is self-consistent. Following the 
discussion after Eq. (3.13), we will estimate Pmax — or Pmax7»7 — |- 

The total viscosity 77 is then the classical plus the new contribution. 



which leads to the an absolute lower bound on the total viscosity: 

/ X53 \ 
V > + ^)W7.j ■ (4.3) 

In a theory with many colors and a weak coupling so that parametrically [40-42] 

e^NlT\ s-^NlT^, rj N^a'^T^ , N^a'^T'^ (4.4) 

the new contribution is parametrically 

W ~ ~ (y'T' ~ (y'N-'r) (4.5) 

Tj 

showing that the new contribution is safely subdominant in any theory which is either weakly 
coupled or has many fields. 

However, in the real world QCD has 3 colors and the coupling at scales of current interest 
is not small. How large is this new contribution to rj in this case? To address this, we need 
estimates for rjci, r^, and e+P = sT for real-world QCD. The easiest is sT, which can be 
determined on the lattice. According to Borsanyi et al [43] Fig. 12, between T — 200MeY 
and T = SOOMeV s/T^ rises from about 10 to about 14. Bazavov et al [44] Fig. 8 gives 
comparable numbers, with s/T^ rising from 11 to 16 in the same range, naturally scales 
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as rjci/sT = 7^, and estimates in QCD vary from 2.67^, the value in strongly coupled Af—A 
SYM [32], to~ 57^, the value in weakly coupled QCD [42]. The trend is that, while r^^/^ is 
a pure number of order unity, it is smaller in more strongly-coupled contexts and larger in 
weakly-coupled contexts. So in the context of QCD with a small shear viscosity it probably 
makes sense to assume Tjt/'^r) is at the low end of this estimated range. Regarding r/ci, the 
lowest estimates for t^ci/s are around 0.08, while r^ci/s = 0.16 may be on the high side in 
terms of fitting elliptic fiow data [10]. 

If we estimate that r^r = 3rici/sT, rid/ s = 0.08, and s = lOT^, we find pmax ~ 2T and 
r^new ~ 0.36T3 - .036s. The true value of rj/s would then be 0.08 + 0.036 = 0.116. Varying 
the value of rjd/s while holding the other estimates fixed, this is close to the minimum value 
of the total viscosity. On the other hand, if rjci — .16s but the other estimates are the 
same, then pmax ~ and 77new ~ .09T^ ~ .01s, which is a negligible correction to the total 
viscosity. Therefore the importance of the "new" contribution to viscosity is quite sensitive 
to the value of rjci/s; for real- world values of other parameters rj/s — 0.08 appears to be 
impossible, but 77/s = 0.16 is not. 

B. Relaxation time 

The presence of a nonanalytic term in the frequency expansion for G^'^^{u) implies that 
hydrodynamics at second order in gradients does not, strictly speaking, work. However, in 
practice we are usually interested in applying hydrodynamics over some range of time, with 
some limited time resolution, and with a finite accuracy tolerance. If the nonanalytic term 
is sufficiently small compared to the term for the frequency scales which are of actual 
importance in a particular problem, then there may be no issue in practice with using second 
order hydrodynamics. To see whether this is the case, we should estimate the frequency c^min 
where the cu^ term is larger than the tIciTt^uj'^ term. Clearly, for u < cUmin, the second-order 
treatment becomes invalid. 

Again using the parametric estimates of Eq. (4.4), the frequency expansion of G^^^^, 
Eq. (3.1) plus Eq. (3.13), is parametrically of form 

~ N^T^co^ - iN^a-^T'u' + (1 + i)N^a'T-2uj'l + N^a-^T^u^ . (4.6) 

The frequency scale where the lu^ term first dominates the cu^ term is parametrically 

7V,VTia;i ~ 7V,2a-^rV =^ u^,^^ N'^^T . (4.7) 

Therefore, in a theory which cither has a large number of colors or has weak coupling, 
the frequency scale where the new term becomes important is parametrically tiny. Then 
there is no obstacle to using second-order hydrodynamics for uj > cUmin; the second-order 
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e+P 




'/d/-s 


^iiiin 


lOT^ 


322ci 
sT 


0.08 


7T 


lOT^ 




0.08 


2.6T 


16T4 




0.08 


2.8T 




•iVcl 


0.16 


T/18 






0.16 


r/125 



TABLE I: Frequency scale Wmin where the uj'^ term equals the rjcYT-j^uP term in G^^^{ui). For 
cj < Wmin, sccond-order hydrodynamics is inconsistent. 



treatment will only become invalid at a frequency scale where it is in any case almost 
irrelevant compared to the uj^ viscous term. 

Now we turn to real- world QCD. To determine the breakdown scale of second-order 
hydrodynamics we again compute the scale where the real part of the term in Eq. (3.13) 
equals the tJciTt^uj'^ term in Eq. (3.1): 



T3 

We illustrate the consequences for real- world QCD in Table I, using some of the estimates 
for s/T^, Ttt, and rjci discussed in the last subsection to evaluate the frequency scale cUmin- 
Below this frequency second-order hydrodynamics is certainly not consistent. What we 
see is that, for reasonable values of s/T^ and of r^, this frequency scale is very low if 
Vci/s — 0.16, but it is very high for rjd/s — 0.08. Therefore, second-order hydrodynamics 
can be applied to QCD if rjci/s ~ 0.16 above the transition temperature; it breaks down 
below frequency scales of order T/20, which corresponds to time scales above 20 Fermi, 
safely above the microphysical scale it is intended to model in heavy-ion coUisions (5-10 
fm). But if 77ci/s ~ 0.08, there is a problem. Second-order hydrodynamics is then already 
inconsistent on frequency scales of order 2.5T, corresponding to time scales above 0.4 Fermi. 
This is a shorter time scale than our estimate p^ax ~ 1/2T for the shortest scale on which 
hydrodynamics is to be reliable. So in this case there would be no range of scales on which 
second-order hydrodynamics is applicable. 



C. Viscosity — again 

Let us accentuate the conclusions drawn in the previous subsection by again considering 
the correction to the viscosity in (3.13). As discussed above, Eq. (3.13) contains a constant, 
Pmax dependent correction to the viscosity, but its actual value depends on the estimate for 

Pmax- 

A different way to bound what values of viscosity are compatible with second-order 
hydrodynamics can be obtained by considering the imaginary part for the co^^^ term, which 
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FIG. 2: Examples for the viscosity over entropy density bound (4.3) (left) and applicability of 
second-order viscous hydrodynamics (4.9) (right). The unknown parameters and Pmax were 
assumed to be Tj^/jn = 3 and Pmax = 1/(2t^) and s/T^ was evaluated using lattice QCD equations 
of state from the hotQCD [44] and Wuppertal-Budapest [43] collaborations. 



can be interpreted as a frequency-dependent correction to the viscosity, 

r]{uj) 



(2407r)7^/^ 

Note that the first and second term become of the same order at u = u„ 

(2407r)2 



where 



(7+(3/2)^ 



Second-order hydrodynamics is certainly no longer applicable at w > Wmax, because the 
frequency dependent contribution to the viscosity becomes large. On the other hand, we 
found above that second-order hydrodynamics also breaks down for uj < Wmin, because there 
the frequency dependent contribution to the relaxation time becomes dominant. 

So for a second-order hydrodynamic description to be applicable at all, there has to be a 
frequency window cjmin < a; < cUmax which can only exist if 



giT(7+(3/2)i)i 
(2407r)i(rj7^)i 



(4.9) 



This is a lower bound on viscosity if second-order hydrodynamics is to have any range of 
validity. If we require the range of allowed frequencies u to be non-zero, the viscosity will 
have be even higher than this bound. 

The temperature dependence of the minimal viscosity can be evaluated when using lattice 
QCD data [43, 44] for the entropy s and using r.n-/7»? = 3, and is shown in Fig. 2. We find 
that the minimal viscosity as defined in Eq. (4.9) varies from about i]/s = 0.09 at 200MeV 
to about Tj/s = 0.07 in the Stefan-Boltzmann limit. These values are comparable to what 
we found is the lowest possible value oir)/ s in Subsection IV A. 
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V. CONCLUSIONS 



We have shown that, just as in non-relativistic hydrodynamics, in a relativistic setting 
it is not self-consistent to consider the hydrodynamic gradient expansion to second order. 
The effects of thermal fluctuations in the hydrodynamic variables themselves contribute to 
the hydrodynamic evolution of the longest wavelength modes. These effects are suppressed 
both in weakly coupled theories and in theories with many degrees of freedom such as QCD 
with many colors; but they can be important in real- world QCD. 

The correction to the shear viscosity is positive and finite, so shear viscosity is a well 
defined quantity in 3+1 dimensions. However, in real- world QCD, if r]c\/s is very small then 
the hydrodynamic fluctuations can make a significant additional contribution; estimating 
their size from our calculation, it does not appear to be possible for the shear viscosity to 
entropy ratio of real-world QCD at T = 200MeV to be smaller than 77/s ~ 0.1. 

The issue becomes more severe for the second-order coefficient r^. Strictly speaking this 
coefficient cannot be defined; its definition assumes an analytic structure for the stress-stress 
correlation function which is violated due to long-wavelength hydrodynamic fiuctuations. In 
practice this is only really an issue at long time scales (low frequencies), where the inclusion 
of second-order effects may not be important anyway. However, how long "long" is depends 
on the value of the viscosity and the entropy density. For real-world QCD, if 7]^/ s ~ 0.16 
then there is a wide range of frequencies where the second-order theory is applicable; it only 
fails at such low frequencies that the difference between first-order and second-order hydro 
is insignificant at the affected scales. In this case there would be no real problem in practice 
with using the second-order theory to e.g. model heavy-ion collisions. But if rjci/s ~ 0.08, 
then there is no range of scales where the application of second-order hydrodynamics is 
consistent in real- world QCD. 

While our calculation focussed on the case of relativistic hydrodynamics, a similar cal- 
culation should be possible in the non-relativistic setting, presumably with qualitatively 
similar results. We leave this interesting project for future work. 

The effects we describe should not be present in dissipative hydrodynamic simulations of 
heavy ion collisions as they are currently conducted. That is because, currently, such simu- 
lations include dissipative viscous effects, but they do not include the fluctuations required 
by the fluctuation-dissipation theorem to ensure that hydrodynamic modes equilibrate with 
mean thermal excitation amplitudes. So the short-wavelength hydrodynamic waves respon- 
sible for the effects we are discussing get quenched in existing hydrodynamic simulations. It 
would be very interesting to try to include thermal fluctuations in hydrodynamic variables 
consistently in hydrodynamic studies of heavy ion collisions. 
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